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Abstract 

We studied previously the impact on light-element abundances of gravitinos decaying during or after 
Big-Bang nucleosynthesis (BBN). We found regions of the gravitino mass 7713/2 and abundance £3/2 
plane where its decays could reconcile the calculated abundance of 7 Li with observation without 
perturbing the other light-element abundances unacceptably. Here we revisit this issue in light 
of LHC measurements of the Higgs mass and constraints on supersymmetric model parameters, 
as well as updates in the astrophysical measurements of light-element abundances. In addition 
to the constrained minimal supersymmetric extension of the Standard Model with universal soft 
supersymmetry-breaking masses at the GUT scale (the CMSSM) studied previously, we also study 
models with universality imposed below the GUT scale and models with non-universal Higgs masses 
(NUHM1). We calculate the total likelihood function for the light-element abundances, taking into 
account the observational uncertainties. We find that gravitino decays provide a robust solution to 
the cosmological 7 Li problem along strips in the (7773/2^3/2) plane along which the abundances of 
deuterium, 4 He and 7 Li may be fit with Xmin ~ 3, compared with \ 2 ~ 34 if the effects of gravitino 
decays are unimportant. The minimum of the likelihood function is reduced to x 2 < 2 when the 
uncertainty on D/H is relaxed and < 1 when the lithium abundance is taken from globular cluster 
data. 
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1 Introduction 



The good agreement of the observed astrophysical abundances of the light elements with 
those calculated in standard Big-Bang nucleosynthesis (SBBN) [1] - jl] imposes interesting 
constraints on late-decaying massive particles, which are generic features of plausible exten- 
sions of the Standard Model, such as supersymmetry. There is even the possibility that the 
decays of massive particles could improve the agreement between calculation and observa- 
tion, particularly in the case of 7 Li [5] - [15] . whose abundance is calculated in standard BBN 
to be considerably higher than that observed in field halo stars or even globular clusters [I] . 



In a previous paper [11] . we studied this possibility in the context of simple supersymmetric 
models with late-decaying massive gravitinos, taking into account the uncertainties in the 
astrophysical observations 0. We found that, for suitable values of the gravitino mass m 3 / 2 
and abundance C3/2 = 777 3 / 2 n 3 / 2 /n 7 , gravitino decays could reduce the combined \ 2 likelihood 
function for the deuterium, 4 He and 7 Li abundances from the SBBN value of ~ 32 to ~ 5.5, 
thus offering a possible solution to the cosmological 7 Li problem. 

Our previous analysis was in the context of the minimal supersymmetric extension of 
the Standard Model with soft supersymmetry-breaking mass parameters constrained to be 
universal at the GUT scale (the CMSSM). Moreover, we chose specific values of the CMSSM 
parameters that were allowed, even favoured, before the advent of the LHC searches for 
supersymmetry. These early parameter choices are now excluded by both the unsuccessful 
searches for supersymmetry at the LHC [16J and the determination of the Higgs mass at 
125-126 GeV [17]. Partly for this reason, there is increased interest in alternatives to the 
CMSSM. 

In this paper, we revisit our previous study including the WMAP9 estimate of the baryon 
density [ISj and, where appropriate, updated estimates of the observational and nuclear- 
reaction uncertainties. In the absence of effects due to gravitino decays, we find that the 
SBBN value of the combined likelihood function \ 2 — 33.7. Here we study the effects of grav- 
itino decays in CMSSM models that are favoured in a recent global analysis of the CMSSM 
parameter space [TO], and extend our study to include models with non- universal Higgs 
masses [2"0"tl2~T] (NUHM1) and models with universality imposed below the GUT scale [2~Tjl2"2"] 
(subGUT). 

In all cases, we find strips in the (7713/2^3/2) plane where the combined \ 2 likelihood 
function may be reduced to ^ 3 if the halo-star estimate of the 7 Li abundance is used, or 
< 1 if the globular-cluster estimate is used (to be compared with an SBBN value of ~ 23.8 
if globular cluster 7 Li is used). If we relax the uncertainty in the D/H abundance (as may 
be warranted by the large dispersion in the data), we find that \ 2 ~ 1-2- We conclude that 
the late decays of massive gravitinos offer a robust solution to the cosmological 7 Li problem. 
The likelihood function for the light-element abundances is minimized for a gravitino mass 
777.3/2 between ~ 4.6 and ~ 6.2 TeV and an abundance £ 3 / 2 in a narrow range between ~ 1.0 
and ~ 2.6 x 10~ 10 GeV in the cases we study. 

■""We also took into account the uncertainties in the relevant nuclear reaction rates, finding them to be 
subdominant, which is also the case here. 
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2 The Observed Light- Element Abundances and their 
Uncertainties 



The ranges we adopt for the observed light-element abundances, follow largely those used in 
our previous work (23]. In that paper we delineated ranges that we considered acceptable, 
problematic, excluded and strongly excluded, as specified in Table [fl which is taken directly 
from [23] . In the subsequent figures, these ranges are left unshaded and coloured yellow, red 
and magenta, respectively. We now comment on the ranges of the light-element abundances 
that we consider to be preferred and use in the subsequent x 2 analysis: more discussion can 
be found in [23] . 
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< 0.24 




< 0.5 




Acceptable 
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[0.3, 1.0] 


[0.24, 0.27] 


< 0.05 


[0.5, 2.75] 


< 0.3 


Problematic 


> 3.7 


> 1.0 


> 0.27 


> 0.05 


> 2.75 


> 0.3 


Excluded 


> 5.0 


> 3.0 


> 0.28 


> 0.1 


> 10 


> 1.0 


Strongly excluded 


> 10 


> 5.0 


> 0.29 


> 0.2 


> 30 


> 3.0 



Table 1: The ranges of light-element abundances whose comparisons with observation we 
consider, following [23], to be acceptable, problematic and (strongly) excluded, as shown in 
the unshaded/yellow/red/magenta regions in the Figures. 



D/H 

We assume the central value [24] 

J = (3.01 ±0.21) x 10~ 5 , (1) 
/ p 

corresponding to the deuterium abundance measured in 11 quasar absorption systems [25] 
and the error in the weighted mean. The SBBN value for D/H at the current baryon-to- 
photon ratio is (2.51 ±0.17) x 10~ 5 . The combined uncertainty (theory and observation), is 
±0.27 x 10~ 5 . We also consider the sample variance in the data, which leads to a combined 
uncertainty of ±0.70 x 10~ 5 that is considerably less restrictive. Using the smaller error 
corresponds to being more conservative in claiming a solution to the 7 Li problem, whereas 
the more relaxed version of the D/H constraint corresponds to a more conservative approach 
to excluding scenarios. We note that we have taken a conservative approach in that we have 
not included the recent measurement of D/H in [26], as the spectrum suffers from the lack of 
a visible Lyman a emission feature since the absorber and emitter are at the same redshift. 



2 



Inclusion of this point would bring the mean downward (to 2.65 x 10~ 5 ) with (in our view) 
an unrealistically small uncertainty of 0.11 x 10~ 5 (and sample variance of 0.36 x 10 -5 ). 
Furthermore, until several such measurements confirm the precision of this value, we cannot 
be assured that no in situ destruction of D/H has taken place, given the dispersion in the 
extant data [24J. 

3 He/D 

As discussed in |23j, we expect the ratio 3 He/D [27] to be a monotonically increasing function 
of time, and therefore consider the solar ratio 3 He/D ~ 1 [28] to be a conservative upper 
bound on the BBN ratio. As we see later, this constraint is satisfied comfortably within all 
the region of interest, so we do not include it in our global x 2 function, and hence do not 
need to assign an uncertainty. 



Following [23], we use the estimate [29] 

Y p = 0.2534 ± 0.0083 (2) 

based on an analysis of the regression of Y vs. O/H. This should be compared to the SBBN 
value of 0.2487 ± 0.0002. 



6 Li/ 7 Li 

As in [23], we assign an upper limit 6 Li/ 7 Li < 0.05, and (like the 3 He/D constraint) do not 
include it in our global \ 2 function. 

7 Li/H 

The cosmological 7 Li problem [3j refers to the discrepancy between the SBBN prediction of 
7 Li/H = (5.1li$; 62 ) x 10~ 10 and observations [HD]- We adopt here the observational range [31] 

C^f) = (1-231JS) x HT 10 , (3) 

V n / halo* 

leading to a combined uncertainty that we take as ±0.71 x 10~ 10 . We note, however, that 
the 7 Li abundance in globular cluster stars is reported to be a factor ~ 2 higher [32l[33] . In 
our analysis below, we also derive results assuming [32] 7 Li/H = (2.34 ± 0.05) x 10~ 10 with 
a combined uncertainty of ±0.62 x 10~ 10 . 

9 Be/H 

We note finally that 9 Be is also observed in halo dwarf stars, with an abundance that varies 
strongly with metallicity, see, e.g., [31]. The observation with the lowest metallicity has 
an [O/H] value of about -2.5 with a 9 Be/H abundance of 3 x 10 -14 . We consider this 
to be a conservative upper limit on the 9 Be/H abundance, while noting that there is one 
observation |35j of 9 Be with an abundance about 3 times lower. As in the cases of 3 He/D 
and 6 Li/ 7 Li, we do not include this observable in our global x 2 function. 
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Finally, we note that our SBBN results are based on the WMAP9 determination of the 
baryon density [18] that corresponds to a baryon-to-photon ratio of 77 = (6.20 ±0.14) x 10~ 10 . 
WMAP9 also provides a value of the cold dark matter density that can be used as an upper 
limit on the possible cosmological density of supersymmetric dark matter neutralino particles 
X produced through thermal freeze out and by gravitino decays: Q x h 2 < 0.1138 ±0.0045. In 
all of the models considered below, the relic density of neutralinos from thermal annihilations 
has been adjusted to yield the central WMAP value of Q x h 2 = 0.1138 using the SSARD 
code [36], so all the neutralinos produced by gravitino decays contribute to the global x 2 we 
define below. 

3 Supersymmetric Models Studied 

We have studied the light-element abundances in ten different supersymmetric models that 
respect the LHC and other current constraints. 

Four of these models are taken from [19] @, which studied the parameter spaces of the 
CMSSM and the NUHM1, taking into account not only the direct searches for supersym- 
metry via missing-transverse-energy events at the LHC [16], but also the mass measured 
for the Higgs boson discovered by ATLAS and CMS [17] and searches for the rare de- 
cay B s — > /i + yU~ [37] H and the XENON100 direct search for dark matter scattering |40j . 
Ref. [19] found two well-defined and almost equally good local minima of the global like- 
lihood function for the CMSSM, with quite distinct values of the scalar and gaugino soft 
supersymmetry-breaking parameters, m and mi/ 2 , and the ratio of Higgs v.e.v.s, tan/3. 
Their values at the local best fits are shown in the first few columns of the first two rows in 
Table |2j The first of these points lies on a strip where the density of the lightest supersym- 
metric particle (LSP), the neutralino Xi is brought into the WMAP range by coannihilations 
with the nearly-degenerate next-to-lightest sparticle, the lighter stau f\. The second point 
is located in a rapid-annihilation funnel where s-channel annihilations through heavy Hig- 
gses control the relic density. We also list the corresponding best-fit values for the trilinear 
soft supersymmetry-breaking parameter, A , at the GUT scale, and note the positive sign 
assumed in [TH] for the Higgs mixing parameter, /z, whose numerical value is fixed by the 
electroweak vacuum conditions in the CMSSM. The next two rows of Table |2] list the cor- 
responding parameters for the best fit adapted from [13] in the NUHM1 framework, and a 
representative good fit with larger values of m ,mi/ 2 and tan/3. In these cases the corre- 
sponding values of fi are not fixed by the electroweak vacuum conditions, and we note their 
fit values. 



2 Their parameters are adapted for study using the SSARD code [55] . 

3 The results in [TO] are very little affected by incorporation of the recent measurement of BR(_B S — > 
by the LHCb Collaboration [38], whose central value is quite close to the combination of pre- 
vious data, and very consistent with the predictions of the CMSSM and NUHM1 fits in [TOJ: see 
http://mastercode.web.cern.ch/mastercode/ for this and other current results from the Mastercode 
group. The best-fit models of [TH] are also compatible with the preliminary constraints from LHC 8-TeV 
data [5S]. 
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ID 
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mi/2 


m 


A 


tan/3 


/i 


m x 


m h 


^3/2 


C3/2 


T3/2 


V 2 . 

A mm 


1 


CMSSM 




US] 


905 


361 


1800 


16 


> 


395 


123.8 


4560 


1.5 x 10" 


-10 


208 


2.81 


2 


CMSSM 




[T9l 


1895 


1200 


1200 


50 


> 


857 


123.3 


5520 


1.8 x 10~ 


-10 


231 


2.86 


3 


NUHM1 




[191 


970 


345 


2600 


15 


2600 


427 


123.8 


4600 


1.2 x 10~ 


-10 


220 


2.82 


4 


NUHM1 




IT91 


2800 


1040 


2100 


39 


3800 


1288 


124.0 


6200 


2.6 x 10" 


-10 


276 


3.14 


5 


CMSSM 


Fie 


2d of I2T1 


1115 


1000 


2500 


40 


> 


496 


124.8 


4800 


1.6 x 10~ 


-10 


213 


2.87 


6 


NUHM1 


Fie 


5b of I2TI 


1175 


1500 


3000 


40 


500 


499 


125.9 


5000 


2.6 x 10" 


-10 


188 


2.86 


7 


NUHM1 


Fie 


6b of I2T1 


1300 


1000 


2500 


30 


-550 


550 


125.5 


4700 


1.0 x 10" 


-10 


258 


2.87 


8 


subGUT CMSSM 


Fie 


9c of [2T1 


2040 


2200 


5500 


10 


> 


1554 


126.7 


5400 


1.6 x 10~ 


-10 


214 


2.96 


9 


subGUT mSUGRA 


Fie. 


10d of [21] 


2400 


4000 


Polonyi 


36 


> 


1099 


125.4 


6000 


1.6 x 10" 


-10 


239 


2.91 


10 


subGUT mSUGRA 


Fig. 


lOd of [21] 


1700 


2000 


Polonyi 


33 


> 


1110 


124.0 


5100 


1.6 x 10" 


-10 


219 


2.89 


11 


CMSSM^ 




m 


905 


361 


1800 


16 


> 


395 


123.8 


4440 


1.5 x 10- 


-10 


230 


1.25 


12 


CMSSM( b ) 




m 


905 


361 


1800 


16 


> 


395 


123.8 


4520 


1.0 x 10- 


-10 


215 


0.52 


13 


CMSSM^ 




m 


905 


361 


1800 


16 


> 


395 


123.8 


4360 


7.1 x 10- 


-11 


245 


0.37 



Table 2: The models studied, with references, their input parameters and the corresponding values of m h calculated using 
FeynHiggs fJTjJ (which have an estimated theoretical uncertainty ^ 1.5 GeV), the best-fit gravitino mass m 3 / 2 and abundance 
(3/2, and the minimum x 2 in a global fit to the observed values of the light-element abundances. All mass parameters 
are expressed in GeV units, and the best-fit lifetime r 3 / 2 is in seconds. The subGUT model in Fig. 9c of WT§ assumes 
M in = 10 9 GeV, and those in Fig. lOd of assume Mi n = 10 10 GeV. In models 11 through 13, the \ 2 and best-fit values 
are computed using (a) the D/H sample variance uncertainty, (b) the 7 Li/H abundance as determined from globular clusters, 
and (c) both the D/H sample variance uncertainty and the globular cluster 7 Li/H. 



We also show in Table [2] the corresponding masses of the lightest neutralino, x-, which is 
assumed in this paper to be the LSP and to constitute (the bulk of) the dark matter. We note 
that m x varies over a large range in the model fits. We also note the values calculated for the 
mass of the Higgs boson, rrih- Considering the theoretical uncertainties in the calculation of 
mh [U], these values are consistent with the LHC measurement rrih ~ 125 to 126 GeV [17] . 

Although the first four rows in Table [2] include the best fits in the CMSSM and NUHM1, 
they do not exhaust the interesting possibilities offered by these models. Accordingly, we 
have also studied three other choices of CMMSM and NUHM1 parameters that are also 
compatible with the LHC constraints and offer complementary possibilities. The fifth row 
lists the parameters of a point towards the tip of the x~^i coannihilation strip for tan/3 = 40 
shown in Fig. 2d of [21], which has larger values of mo, 777,1/2, A) and m x than the best-fit 
point in the first row of Table [2j The sixth row illustrates a new possibility in the NUHM1 
compared to the CMSSM shown in Fig. 5b of [21], namely the appearance at large 777.1/2 of 
a band where the relic LSP density is brought down into the cosmological range because 
the make-up of the x is i n transition between a gaugino-like state at lower 7711/2 which has 
too high a relic density, and a Higgsino-like state at higher mi/2, which has too low a relic 
density. The seventh row in Table [2] is an NUHM1 example with similar values of 777.0,7711/2 
and tan j3, but with ix < 0, which also lies along a transition strip, as seen in Fig. 6b of [21 J . 

Finally, the next three rows in Table [2] illustrate possibilities if the soft supersymmetry- 
breaking parameters are universal at some scale Mj n < Mqut- The eighth row is a CMSSM 
example with M« n = 10 9 GeV, which is located along one of the rapid-annihilation funnels 
visible in Fig. 9c of [21]. The next two are mSUGRA-inspired examples with M in = 10 10 GeV 
that are located on either side of the prominent rapid-annihilation funnel visible in Fig. lOd 
of [21]. In these cases supersymmetry is assumed to be broken by the simplest Polonyi 
superpotential, which specifies the corresponding value of A . These examples populate 
regions of the (777.1/2, T77 ) plane that are not allowed in the CMSSM or the NUHM1 examples 
given above: they all have mi/2 < ttio and large values of m x . 

Since these examples exhibit several features not present in the specific CMSSM mod- 
els studied previously, they serve to probe the robustness of the scenario for resolving the 
cosmological 7 Li problem proposed in [TU[[TT]. 

4 Results 

In this Section we present in some detail results for the LHC-favoured CMSSM and NUHM1 
models presented in the first four rows of Table EJ and then summarize results for the other 
six models. 

In Fig. [1] we display the effects on the light-element abundances of the late decays of a 
gravitino with mass 777 3 / 2 G [2,8] TeV and abundance C3/2 £ [10~ 14 , 10~ 7 ] GeV. We show 
blocks of six panels for each of the models described in the first four rows of Table [2j lower- 
mass CMSSM (upper left), higher- mass CMSSM (upper right), lower- mass NUHM1 (lower 
left) and higher- mass NUHM1 (lower right). For each model, the upper left panel displays 
the effect of gravitino decays on the D/H ratio, the upper middle panel the effect on the 
3 He/D ratio, the upper right panel the effect on the 4 He abundance, the lower left the effect 
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on the 6 Li/ 7 Li ratio, the lower middle the effect on the 7 Li/H ratio, and the lower right the 
effect on the 9 Be/H ratio. In each panel, the unshaded, yellow, red and magenta regions 
correspond to the ranges specified in Table [TJ Solid shadings are used for regions with excess 
abundances, and hashed shadings for regions with low abundances. 

The corresponding panels for the different models exhibit several common features. We 
note first that the D/H constraint typically excludes a large triangular region at small 777.3/2 
and large £3/2; but has no impact in the triangular region at large 7773/2 and small £3/2- Next 
we note that the 3 He/D constraint has no impact in any of these models. The impact of 
the 4 He constraint is always limited to a small region at relatively small 7773/2 and large £3/2 
that lies within the region already excluded by the D/H ratio. The 6 Li/ 7 Li ratio typically 
excludes a band at smaller 777,3/2 that is cut off at very small £3/2- The 7 Li/H ratio excludes 
not only a region at small 7773/2 and large £3/2, but also a region at large 7773/2 and small 
£3/2. In between, there is a banana-shaped region extending from 7773/2 ~ 3 TeV upwards, 
where the 7 Li abundance falls within acceptable limits. The yellow shading at lower right in 
each of these panels reflects the existence of the cosmological 7 Li problem in standard BBN, 
and the unshaded banana indicates that this may be solved by the late decays of massive 
gravitinos. Finally, we note that the 9 Be/H constraint has impact only at small 777.3/2 and 
large (3/2, in a region already excluded by the D/H ratio. 

It is clear from this survey that, in all these four models, the most important constraints 
are those from the D/H and 7 Li/H ratios, and the issue is whether, within the unshaded Li 
banana, there are model parameters where these constraints can be reconciled adequately. 

This issue is addressed in Fig. [2] for the CMSSM and NUHM1 models described in the 
first four rows of Table [2j As in Fig. [TJ the upper left panel is for the lower-mass CMSSM 
model, the upper right panel for the higher-mass CMSSM model, the lower left panel for the 
lower-mass NUHM1 model, and the lower right panel for the higher-mass NUHM1 model. 
In each panel, we display the values of D/H and 7 Li/H for a grid of points in the (777.3/2, C3/2) 
plane. The pattern of points in Fig. [2] arises as follows. Each of the loops (best visible at 
moderate D/H and 7 Li/H) drops nearly vertically, goes through a minimum, and then rises 
towards high 7 Li/H and low D/H. A single loop represents a single value of £3/2- As seen 
in Fig. [TJ scanning from low to high 777.3/2 at fixed £3/2 generally takes us from high D/H to 
low D/H, and from high 7 Li/H to low 7 Li/H to somewhat higher 7 Li/H. This leads to the 
loop behavior. For example the top "row" of blue points seen in the two panels on the left, 
correspond to 777,3/2 = 2 TeV and £ 3 / 2 increases from left to right. We see that the points 
examined in this scan grazes the inner ellipse representing the combined one-cr excursion of 
the observed D/H and 7 Li/H ratios. 

The predominant feature in Fig.[2]is the correlation between lower 7 Li/H and higher D/H 
that is a common feature of many potential solutions to the 7 Li problem [24|. This is in 
fact a welcome feature as some enhancement in D/H is warranted when comparing the mean 
value from observations and the SBBN result. The red crosses represent the best-fit points 
in each of the four models, and lie close to the one-cr ellipses in all cases. We learn from 
these plots that the D/H and 7 Li/H observations can indeed be reconciled, with the D/H 
and 7 Li/H ratios lying within their combined one-tx observational band, somewhat above the 
central values given in Section 2. 
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Figure 1: Plots of the effects on the light-element abundances of the decays of a gravitino with mass m 3 / 2 £ [2,8] TeV 
and abundance ( 3 / 2 G [10~ 14 , 10~ 7 ] GeV. The results are displayed in blocks of six panels for each of the models described 
in the first four rows of Table In each block, the upper left panel displays the effect on the D/H ratio, the upper middle 
panel that on the 3 He/D ratio, the upper right panel that on the 4 He abundance, the lower left that on the 6 Li/ 7 Li ratio, the 
lower middle that on the 7 Li/H ratio, and the lower right that on the 9 Be/H ratio. The unshaded regions in the panels are 
those allowed at face value by the ranges of the light- element abundances reviewed in Section 2, whilst the yellow, red and 
magenta regions correspond to progressively larger deviations from the central values of the abundances, as summarized in 
Table\J\ Solid shadings are used for regions with excess abundances, and hashed shadings for regions with low abundances. 



CMSSM 
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CMSSM 
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Figure 2: Planes giving overviews of the D/H and r Li/H abundances obtainable in the models 
described in the first four rows of Tabled The blue points show the values found in scans of 
the models for different values of m 3 / 2 and ( 3 / 2 . The ellipses represent the one- and two-cr 
regions found by combining the D and 7 Li constraints reviewed in Section 2. The red cross 
marks the best fit found in each model. 

The broader issue of compatibility with all the constraints from light-element abundances 
is addressed in Fig. [3j which explores the mutual compatibility of all the light-element 
constraints in the (TO3/2, C3/2) planes for the models described in the first four rows of Table|2j 
In each case, the magenta region is where one or more calculated abundance lies within a 
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'strongly-excluded' range specified in Table [TJ whereas in the red region the largest deviation 
is in an 'excluded' range, and in the yellow region the largest deviation is in a 'problematic' 
range. In each panel, we see that there is a strongly excluded magenta region at smaller 
777.3/2 and larger (3/2 (where there are deviations for D/H and some other abundances) and 
a disfavoured region at larger 717,3/2 and smaller £ 3 / 2 (where the 7 Li abundance is similar to 
that in standard BBN) separated by an unshaded banana-shaped strip where no calculated 
light-element abundance disagrees significantly with observation. While the exact location 
and shape of the unshaded banana differs in each model, its presence is a generic feature we 
see not only in the four models displayed but in all ten we have studied. 

As a refinement of this analysis of the compatibility between the different light-element 
constraints, we display in Fig. @] contours of the global \ 2 function in the (7773/2, C3/2) planes 
for the models listed in the first four rows of Table [2j The value of \ 2 is determined from 

2 / / y p - 0.2534 \ 2 /D/H — 3.01 x 10 5 \ 2 / 7 Li/H — 1.23 x 10~ 10 \ 2 f^ /2) h 2 \ 2 
X = V 0.0083 J + V 0.27 x 10- 5 ) + V 0.71 x lO" 10 J + y 0.0045 J ' 

(4) 

where 

^ 3/2) = ^^Cs/ 2 (5) 

is the density of neutralinos produced in gravitino decays, p c being the critical density of 
the Universe Q In each case, we see a narrow valley where \ 2 < 6.0 (blue line), representing 
an acceptable solution to the cosmological 7 Li problem. In each case, the rise to \ 2 = 9-2 
(magenta line) is quite rapid, but the further rise at smaller 7773/2 and larger ( 3 / 2 (due to a 
combination of several constraints) is much more rapid than that at larger 777 3 / 2 and smaller 
(3/2 (where it is due solely to the cosmological 7 Li problem). At high £ 3 / 2 the WMAP 
constraint plays a role in closing the contours for \ 2 = 6 and 9.2 as well as the bend and 
flattening of the \ 2 = 33.7 and 50 contours respectively. The black crosses in Fig. H] indicate 
the best-fit points in the different scenarios: their properties (7773/2, C3/2, T 3/2, Xmin) are given 
in the first four rows of Table [2J We see that the best-fit gravitino mass varies between 
4.6 and 6.2 TeV, and its abundance in the narrow range between 1.0 and 2.6 x 10~ 10 GeV. 
The best-fit gravitino lifetimes fall in an even more narrow range, T3/2 ~ 210 — 280 sec. 
Thus the models all show a close similarity in the best-fit gravitino abundance and lifetime 
values. This behavior is typical of decaying particle models we and other have studied, and 
indicates that these parameters exert the most sensitive control on the abundances of the 
light elements. 

All four of the models have a Xmin °f about 3, which represents a significant reduction from 
the SBBN value of 33.7 0. However, there is also a price to be paid, in that we have added 

4 We note that it would in principle be possible to adapt the parameters of the models listed in Table [2] 

so that the density of thermally-produced neutralinos is below the central WMAP value, in which case the 
(3/2) 

constraint on would be relaxed. However, since this constraint is not important in the neighbourhoods 

of the best-fit points, we have not explored this option. 

5 The differences between these values and those found in [TT] are due to a combination of factors, including 
the updated D/H abundance as well as the revised value of fi^/i 2 as well as the new supersymmetric model 
parameters. 
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Figure 3: Planes exploring the compatibility of all the light- element constraints in the 
{ m 3/2, C3/2) planes for the models described in the first four rows of Tabled In each case, the 
magenta region is where one or more calculated abundance lies within a 'strongly- excluded' 
range specified in Tabled whereas in the red region the largest deviation is in an 'excluded' 
range, and in the yellow region the largest deviation is in a 'problematic' range. In the 
unshaded banana-shaped regions, there is no significant deviation in any of the calculated 
abundances. 
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2 parameters (777.3/2 and C3/2) to achieve that benefit. As a result we have also decreased the 
number of degrees of freedom from 3 to 1. Nevertheless, the fit probability for x 2 per degree 
of freedom of 3/1 is far superior to that for 33.7/3. 

As also seen in Table EJ we have explored the effects of late-decaying massive gravitinos in 
six other scenarios: three other CMSSM and NUHM1 scenarios and three subGUT scenarios, 
all of them consistent with the current LHC and other constraints and representing different 
ways in which the relic x density can be brought into the range favoured by WMAP and 
other experiments, as discussed earlier. The main features of our results for these scenarios 
are similar to those of the previous models shown and we list the parameters of the best 
fits in the 5th to the 10th row of Table [3 Each of them lies within the ranges found in the 
models discussed above. 

We have studied effects in the first of the models in Table [2] under different assumptions 
on the D/H and 7 Li/H abundances, with the results shown in Fig. [51 In the left panel we 
show contours of the x 2 function if the D/H constraint (PQ) is relaxed, keeping the same 
central value but assuming an uncertainty of ± 0.70 x 10~ 5 , given by the sample variance. 
In this case, the banana-shaped region with x 2 < 6.0 is significantly expanded, and the 
minimum value of x 2 — 1-25 (compared with the SBBN value of 30.7). The middle panel of 
Fig. [5] shows results obtained keeping the more restrictive uncertainty in D/H but assuming 
the globular cluster value for 7 Li/H = (2.34 ± 0.62) x 10~ 10 . In this case, the 'banana' is 
even broader, and the minimum value of x 2 — 0.52 (compared with the SBBN value of 23.8). 
Finally, the right panel shows the effects of relaxing both the D/H and 7 Li/H constraints, 
and we see that the combined result is more similar to the result of relaxing the 7 Li/H 
constraint alone than that of relaxing the D/H constraint alone. In this case the minimum 
value of x 2 = 0.37 (compared with 20.8 in SBBN), demonstrating again that relaxing the 
7 Li/H constraint is more important than relaxing the D/H constraint. 

Table [3] presents light-element abundances predicted at the best-fit points for each of 
the models defined in Table [2j We see that all models select very similar best-fit values 
of D/H and 7 Li/H, in line with the trends found in Fig. |2l namely a combination of the 
highest tolerable D/H and lowest tolerable 7 Li/H. We also see that, in all models, there 
is no change in the helium abundance from its SBBN value. Finally, Table |3] shows the 
best-fit prediction for the relic x abundance due to gravitino decays, f^ 3 /i 2 . These vary 
from ~ (0.3 — 2) x 10 -3 , and thus are at most about half of the observational uncertainty 
in f2 C( j m /i 2 . The relic abundance constraint is therefore not a substantial contributor to the 
global x 2 function in the regions of the best-fit points. 

5 Summary and Conclusions 

We had shown previously [11] that the late decays of massive gravitinos offered a possible 
solution of the cosmological 7 Li problem in a number of benchmark CMSSM scenarios pro- 
posed before the start-up of the LHC. The subsequent non-discovery of supersymmetry at the 
LHC [IB], the discovery of a Higgs boson [T7], the measurement of B s — > fi + fi~ [37] and the 
XENON100 dark-matter search experiment [10] have squeezed severely the allowed region 
of the CMSSM parameter space [19], and stimulating the exploration of alternative super- 
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CMSSM tan|3=16 



CMSSM tan|3=50 




Figure 4: Contours of the \ 2 function in the (1713/2, C3/2) planes for the models described in 
the first four rows of Table\B' X 2 — 6.0 (blue), 9.2 (magenta), 33.7 (red) and 50 (black). 
Also marked by black crosses are the best-fit points in the various models. 
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Figure 5: Contours of the x 2 function in the (rn.3/2, C3/2) planes for the low-mass CMSSM 
model described in the first row of Table [H under different assumptions on the D/H and 
7 Li/H abundances. In the left panel we relax the D/H constraint by using an increased 
uncertainty ±0.70 x 10~ 5 ; in the middle panel we use the globular cluster value for 7 Li/H 
= (2.34 ±0.05) x 10~ 10 ; and in the right panel we combine these assumptions. The contours 
are coloured as follows: x 2 = 2.3 (green), x 2 = 6.0 (blue), 9.2 (magenta), and 50 (black). 
The red contour corresponds to the SBBN x 2 for each case, namely 30.7, 23.8, and 20.8 for 
the left, middle, and right panels respectively. Also marked by black crosses are the best-fit 
points under these different assumptions. 

symmetric scenarios such as NUHM1 models with non-universal Higgs masses and subGUT 
models in which soft supersymmetry-breaking masses are assumed to be universal at some 
scale below the GUT scale [2Tj . In this paper, we have explored the most favoured sets of 
CMSSM and NUHM1 parameters not yet excluded by the LHC [19], as well as a number of 
alternative models within the NUHM1 and subGUT frameworks, as summarized in Table [2j 

In contrast to the constraints on supersymmetric model parameters, the constraints on 
primordial light-element abundances have remained essentially unchanged since [11J, as sum- 
marized in Table [TJ Thus the cosmological 7 Li problem remains. 

In this paper we have repeated the global x 2 likelihood analysis of [11] , and studied the 
regions of the gravitino parameters (^3/2^3/2) that offer the best prospects of solving the 
cosmological 7 Li problem. In all the models studied, we find global Xmin ~ 3 in the (m 3 / 2 , C3/2) 
plane and therefore the cosmological 7 Li problem can be regarded even more effectively solved 
than in [II] . Typical ranges of the gravitino parameters are 4.6 TeV < m 3 / 2 < 6.2 TeV and 
1.0 x 10- 10 GeV < C3/2 < 2.6 x lO" 10 GeV. 

In view of this persistence of the massive gravitino solution despite the impact of LHC 
results on the supersymmetric parameter space, we conclude that the late decays of massive 
gravitinos provide a robust solution to the cosmological 7 Li problem. 
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ID 




10 5 D/H 


10 10 7 Li/H 




v 2 . 
z\ xti i n 


1 


24? 


37 


3.27 


2.12 


5.0 x 10- 


-4 


2.81 


2 


D 24? 


]7 


3.28 


2.09 


1.1 X 


io- 


-3 


2.86 


3 


n 24? 


]7 


3.26 


2.14 


4.4 x 1(T 


-4 


2.82 


4 


n 24? 


]7 


3.29 


2.11 


9 1 v 


1U 


-3 


3.14 


5 


n 24? 


j i 


3.32 


2.01 


6.5 x 


io- 


-4 


2.87 


6 


n 94? 


\7 


3.27 


2.11 


1.0 x 


io- 


-3 


2.86 


7 


n 24? 


\7 


3.29 


2.08 


4.7 x 


io- 


-4 


2.87 


8 


n 24? 


\7 


3.25 


2.16 


1.8 x 


io- 


-3 


2.96 


9 


0.24? 


17 


3.31 


2.04 


1.2 x 


io- 


-3 


2.91 


10 


0.24? 


17 


3.28 


2.09 


1.4 x 


io- 


-3 


2.89 


11 


0.24? 


17 


3.55 


1.63 


5.1 x 


io- 


-4 


1.25 


12 


0.24? 


17 


3.10 


2.50 


3.5 x 10- 


-4 


0.52 


13 


0.24? 


17 


3.15 


2.40 


2.5 x 10- 


-4 


0.37 



Table 3: Predictions at the best-fit points for the models defined in Tabled for light- element 
abundances and the neutralino abundance fl^ h 2 arising from gravitino decays. 
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